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A generation of magnetic fields from cosmological density perturbations is investigated. In the 
primordial plasma before cosmological recombination, all of the materials except dark matter in the 
universe exist in the form of photons, electrons, and protons (and a small number of light elements). 
Due to the different scattering nature of photons off electrons and protons, electric currents and 
electric fields are inevitably induced, and thus magnetic fields are generated. We present a detailed 
formalism of the generation of cosmological magnetic fields based on the framework of the well- 
established cosmological perturbation theory following our previous works. We numerically obtain 
the power spectrum of magnetic fields for a wide range of scales, from k ~ 10~ 5 Mpc -1 up to 
k ~ 2 x I0 4 Mpc -1 and provide its analytic interpretation. Implications of these cosmologically 
' generated magnetic fields are discussed. 

PACS numbers: 91.25.Cw, 98.80.-k 

I. INTRODUCTION 

> 

' There are observational evidences which indicate that magnetic fields exist not only galaxies, but also in even larger 
, systems, such as cluster of galaxies and extra-cluster spaces In addition to gravity, they would play an important 
t— I ■ role in the formation processes of various objects and their dynamical evolution in the universe. The effects of large 
scale magnetic fields on galaxy formation and on formation of the first stars are extensively studied 0, [H, 0, d, @| . 

Yet, the origin of such large scale magnetic fields is still a mystery [lj. It is now widely believed that the magnetic 
fields at large scales are amplified from a tiny field and maintained by the hydro-magnetic processes, i.e., the dynamo. 
However, the dynamo needs a seed field to act on and does not explain the origin of magnetic fields. As far as the 
Oh magnetic fields in galaxies are concerned, the seed fields as large as io~ 20 ~~ 30 G are required in order to account for 
the observed fields Q of order \\i G at the present universe. Discoveries of magnetic fields at damped Lya systems 
Q in the early universe and an intervening galaxy at intermediate redshift Q may require even larger seed fields. 

The question arises here is, thus, what is the origin of such seed fields? Most of previous theories to explain 
the origin of seed fields can be categorized into two, i.e., astrophysical and cosmological mechanisms. Almost all 
astrophysical mechanisms of generating seed magnetic fields are based on the Biermann battery effect [Toj . This 
mechanism is operative when the gradients in thermodynamic quantities, such as temperature and density, are not 
parallel to that of pressure. This mechanism has been applied to various astrophysical systems, which include stars 
supernova remnants Il 2| . |l3ll, protogalaxies [14j . large-scale structure formation |15| . and ionization front at 
cosmological recombination fl6| . These studies show that magnetic fields with amplitude 10~ 16 ~ 10~ 21 G could 
be generated. Recently, it is proposed that the Weibel instability at the structure formation shocks can generate 
rather strong magnetic fields 10 -7 G (UIHI- However, the coherence-length of seed fields generated by astrophysical 
mechanisms, especially the Weibel instability, tends to be too small to account for galaxy-scale magnetic fields. 

On the other hand, cosmological mechanisms based on inflation can produce magnetic fields with a large coherence 
length since accelerating expansion during inflation can stretche small-scale fields to scales that can exceed the causal 
horizon. However, in the simplest models with the usual electromagnetic field, which is conformally coupled to gravity, 
the energy density in a fluctuation of the field diminishes as a -4 (where a is cosmic scale factor) to lead a negligible 
amplitude of magnetic fields at the end of inflation. Therefore, to create enough amount of magnetic fields which 
survive an expansion of the universe until today, some exotic couplings between the electromagnetic field to other 
fields, such as dilaton [H, US], Higgs type scalar particles [2l|, or gravity [HI], must be introduced. Although some 
models of them are considered to be natural extensions of standard particle model, the nature of generated magnetic 
fields depends highly on the way of extension, whose validity can never be tested experimentally by terrestrial particle 
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accelerators. Moreover, it is recently argued that almost all models that generate magnetic fields during inflation at 
the galactic scale (A ~ 0.1 Mpc) are severely constrained as B x < 10 -39 G [H, [24|, for a blue primordial spectrum. 
This is because anisotropic stress of magnetic fields would produce a large amount of gravitational waves, which then 
spoil the standard Big Bang nucleosynthesis by bringing an overproduction of helium nuclei (and deuterium). 

In addition to the two categories described above, however, there is a third category for the generation of large-scale 
seed fields: cosmological fluctuations in the universe can create magnetic fields prior to cosmological recombination. 
This category dates back to Harrison (1970) [25| and has been rigorously studied in recent years. The attractive point 
is that mechanisms based on cosmological perturbations are much less ambiguous than the previous two categories 
because cosmological pertrubations have now been well understood both theoretically and observationally through 
cosmic microwave background (CMB) anisotropy and large scale structure of the universe. Thus it is possible to make 
a robust quantitative evaluation of the generated magnetic fields. 

Originally, Harrison found that the vorticity in a primordial plasma can generate magnetic fields. This is because 
electrons and ions would tend to spin at different rates as the universe expands due to the radiation drag on electrons, 
arising rotation-type electric current and thus inducing magnetic fields. It was found that magnetic fields of 10 -21 G 
would result if the vorticity is equivalent to galactic rotation by z = 10. More recently, magnetic field generation at 
recombination was investigated by evaluating the induced electric field which arise in the matter flow dragged by a 
dipole photon field, resulting in a seed field 10 -26 G [26|. Following Harrison's idea, Matarrese reported that 10 
G seed fields can arise from the vorticity generated by second-order density perturbations at recombination [27l |. In 
a similar analysis, a larger value, 10 -21 G, was obtained by considering earlier periods when the energy density of 
photons was larger and the photon mean free path was smaller (28|. Other specific second order effects from the 
coupling between density and velocity fluctuations on the seed field generation was evaluated in [2!| . 

In [30|| , we presented a formalism to calculate the spectrum of seed fields based on the cosmological perturbation 
theory. There are three essential points in order to obtain the correct spectrum. First of all, we have to treat electron, 
proton and photon fluids separately to evaluate the amount of electric current and electric field. Secondly, we need a 
precise evaluation of collision terms between three fluids, which generate the difference in motion between them. At 
this point, we found that not only the velocity difference between electron and photon fluids, which is conventionally 
studied, but also photon anisotropic stress contributes to the difference in motion between electron and proton fluids. 
Finally, the second-order perturbation theory is necessary because vector- type perturbations, such as vorticity and 
magnetic field, are known to be absent at the first order. Based on this formalism, in [31j, we calculated the spectrum 
of generated magnetic fields at a range between k ~ 10 -5 Mpc -1 and k ~ 10 Mpc -1 . Then we showed that the 
magnetic field is contributed mainly from the velocity difference between electrons and photons, which we called the 
" baryon-photon slip term" , on large scales (k < 1 Mpc -1 ), while photon anisotropic stress is dominant on small scalles 
(k > 1 Mpc -1 ). As a whole, the spectrum is small-scale-dominant and the amplitude is about 10 -20 G at k = 1 
Mpc -1 . 

In this paper, we present details of our formalism and extend it in several directions. First, in section II, we describe 
the derivation of equations of motion for electron and proton fluids, carefully evaluating the collision term between 
electrons and photons, and then obtain the evolution equation for magnetic fields. We calculate and interpret the 
resulting spectrum numerically in section III. In section IV, we derive the analytical expression of the spectrum, which 
allows us to confirm the numerical results and extrapolate the spectrum into much smaller scales than can be obtained 
numerically. In section V, we discuss about the magnetic felicity, which in fact explicitly vanishes. Finally, in section 
VI, we give some discussions and conclusions. 



II. BASIC EQUATIONS 

Here we derive basic equations for the generation of magnetic fields, i.e., perturbation equations of photon, proton 
and electron fluids. While protons and electrons are conventionally treated as a single fluid, however, it is necessary 
to deal with proton and electron fluids separately in order to discuss the generation of magnetic fields. Let us begin 
with the Euler equations. Those are given by 

m p nu^u pi .^ - enu^Fi^ = C* e + C* 7 , (1) 
m e nu%u ei . ttt + enu£F iM = C e f + C e ? , (2) 

where TO p ( e ) is the proton (electron) mass, u p r e \ is the bulk velocity of protons (electrons), F^i is the usual Maxwell 
tensor. The thermal pressure of proton and electron fluids are neglected. Here n, v = 0, 1,2,3 and i — 1,2,3. The 
r.h.s. of Eq. (Q]) and ^ represent the collision terms. The first terms in Eqs. ([T]) and © are collision terms for the 
Coulomb scattering between protons and electrons, which is given by 



-(u pi - u ei )e 2 n 2 -q , 



(3) 
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where 



7re 2 me /2 , , nt __ 1fi (\ + z\ ~ 3/2 /In A 



,-^MA~ M xlO-V=^j ^j. (4) 

is the resistivity of the plasma and In A is the Coulomb logarithm. As is well known, this term acts as the diffusion 
term in the evolution equation of magnetic field. The importance of the diffusion effect can be estimated by the 
diffusion scale, 

/ T \ 1/2 

A diff = V^F - 100 (j^t) au > (5) 

above which magnetic field cannot diffuse in the time-scale r. Here Ho = 70km/s/Mpc is the present Hubble 
parameter. Thus, at cosmological scales considered in this paper, this term can be safely neglected. 

The other terms expressed by C p ^ 7 are the collision terms for Compton scattering of protons (electrons) with 
photons. Since photons scatter off electrons preferentially compared with protons by a factor of (m e /m p ) 2 , we can 
safely drop the term C p ^ from the Euler equation of protons. This difference in collision terms between protons and 
electrons ensures that small difference in velocity between protons and electrons, that is, electric current, is indeed 
generated once the Compton scattering becomes effective. 



A. Compton Collision Term 

Let us now evaluate the Compton scattering term. In the limit of completely elastic collisions between photons 
and electrons, this term vanishes. Typically, in the regime of interest in this paper, very little energy is transfered 
between electrons and photons in Compton scatterings. Therefore it is a good approximation to expand the collision 
term systematically in powers of the energy transfer. 

Let us demonstrate this specifically. We consider the collision process 

7(Pi) + e-(« i )->7(Pi) + e-(gO ) ( 6 ) 

where the quantities in the parentheses denote the particle momenta. To calculate this process, we evaluate the 
collision term in the Boltzmann equation of photons: 

xS^[p + q-p> - q'}S[E(p) + E e (q) - E(p') - E e (q')] 

X [f e (q')f(p')(l + f(p)) - fe(q)f(p)(l + f(p'))] , (7) 



where f(p) and f e (q) are the distribution functions of photons and electrons, E e (q) = y/q 2 + m 2 is the energy of 
an electron, and the delta functions enforce the energy and momentum conservations. We have dropped the Pauli 
blocking factor (1 — f e ). The Pauli blocking factor can be always omitted safely in the epoch of interest, because f e 
is very small after electron-positron annihilations. Note that the stimulating factor can be also dropped because this 
does not contribute to the Euler equation. 
Integrating over q', we obtain 

cm] = -f — f d3p ' — i^i 2 

W/WJ pj {2n)*2E e (q)J (2n)*2E(p>) 2 E e (\q+ p- p>\) 1 ' 
xS[E(p) + E e (q) - E(p') - E e (\q+p-p'\)] 

x [f e (q + P- P')f{p') - fe(q)f(p)} ■ (8) 

In the regime of our interest, energy transfer through the Compton scattering is small and can be ignored in the 
first order density perturbations. As we already discussed earlier, however, it is essential to take the second order 
couplings in the Compton scattering term into consideration for generation of magnetic fields. Therefore we expand 
the collision term up to the first order in powers of the energy transfer''', and keep terms up to second order in density 
perturbations. 



[f] However, we shall keep up to the second order terms for the purpose of reference. 
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The expansion parameter is the energy transfer, 



E e (q)-E e (q') 



(q + p- p'f (p' - p) • q (p- p') 



r\2 



2m e 



2m e 



m. 



2m e 



(9) 



over the temperature of the universe. Employing p ~ T, we can estimate the order of this expansion parameter as 
O(^pjt) ~ O(-), which is small when electrons are non-relativistic. Note that, in the cosmological Thomson regime, 

2 

electrons in the thermal bath of photons are non-relativistic, p ~ 2m - ' an< ^ the ener gy °f photons is much smaller 
than the rest mass of a electron, p <C m e . Thus, it also holds that q ~ y/2m e p 3> p, and the second term in Eq.([9]) is 
usually smaller than the first one. 

Now let us divide the collision integral into four parts, i.e., the denominators of the Lorentz volume, the scattering 
amplitude, the delta function and the distribution functions, and expand them due to the expansion parameter defined 
above. First of all, the denominator in the Lorentz invariant volume can be expanded to 



E e (q)E e (\q+p-p'\) 



2m 



-If 



-£, 



1 



2m, 



(10) 



where 



m " mi 



(p-p') ■ q 



£ (-E-)2 



ip-p'f 
2ml 



(11) 



Secondly, we consider the matrix element. The matrix element for Compton scattering in the rest frame of the 
electron is given by, 



\M\ 2 
cos/3 



2(47r) 2 a 2 



P-P , 



P_ 

P p' 



P .23 
= — sin p 



(12) 



where p and p' are the energies of incident and scattered photons, p and p' are the unit vectors of p and p , respectively, 
denoting the directions of the photons in this frame. The Lorentz transformation with electron's velocity (q/m e ) gives 
the following relations, 



p _ yjl - (q/m e ) 2 
P 



P^P 



1 — p ■ q/pm e 
Pa- 



using these relations, we evaluate the matrix element in the CMB frame as [32 1 

\M\ 2 = 2(4?r) 2 a 2 (m +Mx +7W ( _a_)2 + M { jjp.) + M^y^ 



Mo = 1 + cos 2 P , 
Mj_ = -2cos/3(l - cos/3) 



P + P' 



A^(^_)2 = 2cos/3(l - cos/3) — - 

e 772. „ 



Mjul = (l-cos/3)(l-3cos/3) 



X ( j^)2 = (1-cos/3) 2 At 



— -(p + pO 

m P 



2cos/3(l - cos/3) 



(q-p)(q-p') 



itif 



(13) 
(14) 

(15) 
(16) 
(17) 

(18) 

(19) 

(20) 
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Thirdly, we expand the delta function to 



6\p-p' + E e (q)-E e (q')] 

X( ^ 35{p - p' + E e (q) - E e {.q')) , 



+ 2 



dE e (q>) 

ld 2 5(p-p' + E e (q) - E e (q')) . 



5{ P -p') + 
5{ P -p') + 



dp 2 

dS(p — p') 
dp' 

dS(p — p') 
dp' 



e (g)=B e (9') ( E e(q) - E e (q')Y 



(p-p') ■ q (p-p 1 ) 2 



where 



V 



2m e 



(p-p') ■ q 



d 2 5(p - p') 
dp 72 



2~<£> 9p' 2 

(p-p') 2 



2m„ 



Finally, the distribution of the electron can be expanded to 



8 f -. B 2 f 

fe(q+p- P ') « /e(€) + ■ (P-P') + (f ~P' l )d^j(P J ~P )■ 

We assume that the electrons are kept in thermal equilibrium and in the Boltzmann distribution: 



f e (q) = n e 



2tt 

m e T R 



3/2 



exp 



(q - m e Vef 



2m e T K 



(21) 



(22) 



(23) 



(24) 



where v e is the bulk velocity of electrons. The derivatives of the distribution function with respect to the momentum 
are given as 



dfe _, f-^ q- m e v e 

oq rn e l e 

d 2 f e df e (q) q l - m e v i 

dqidq.j dqj m e T e 

By substituting above equations, Eq. ([23|) is written as 



fe(q); 



(25) 
(26) 



fe(q + p-P) « fM - f e (q) * ■ (P-P 1 ) + \fM 

m e l e 2 



(J>-P')(q- m-eVe) 
m e T e . 



\p-pf 

m e T e . 



fM 



If. 



(£) ~^^) 



Therefore, we have 

fe(q+p-p')f(p') - fe(q)f(f) = fe(q) \f(p') - f(p)} - f(p')fe(q)F { ±) - f(p')fe(q) 



(27) 



(28) 



Fortunately, it has been known that the leading term (zeroth order term), obtained by multiplying together the 
first term in the delta function and the zeroth order distribution functions, is zero. It means that we only have to 
keep up to the first order terms when we expand the matrix element and the energies, in order to keep the collision 
term up to the second order [33| . 

Therefore we have, 



IMF 



E e (q)E e (q') 



6tt(Tt(Mo +Mji_) 



(29) 
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Combining altogether, we obtain the collision term expanded with respect to the energy transfer as (note that this 
expansion is not with respect to the density perturbations) 

(ftp 1 (fig 

[(Oth order term) + (1st order terms) + (2nd order terms)] + ■ • • 



where 



°w-;(S) 



Oth order term: 



1st order terms: 



(2tt) V (2tt) 3 



6k(TtM 8(p - p')f e (q) f(p>) - f(p)) 



6ira T M f e (q) 



+67ra T MjLf e (q)S(p-p') fiP 1 ) ~ fip) 



2nd order terms: 



6ira T Mofe{q) 

d5(p-p') 
+ dp' V % 

+Q-KCF T MjLfe{q) 



5(P-P')f(p r ) 

fip') - fip) 



dp' 



fip') - fiP) 



(30) 



(31) 



(32) 



l d 2 5(p-p') v2 
2 dp' 2 * 



fiP') ~ fiP) 



d8(p- P >) 



dp' 



Vjl+Vjl) fip')TjL 



S(p-p')f(p')^ + 



dS(p — p') 
dp' 



Vp. 



fiP' ~ fiP)) 



(33) 



From now on, we omit the second order terms. These terms are not only much smaller than the first order terms but 
also may not contribute to the Euler equation at all (sec 34]). Evaluating the first moment of the above collision 
term, we obtain the Compton scattering term in the Euler equation (|41[) as 

d 3 p 



C 



(2tt) 3 
4aT/0 7 an e 



PiC[f] 



Here moments of the distribution functions are given by 

P Pfjif) = Py, 



(2tt)3^ 7 
d 3 P , 4 



(2tt) 3 
d 3 q 



f e (q) = n e , 

3 qifeiq) = PeVei, 



(27T) 

d 3 p 1 1 

j^r S p- 1 PiPjf~,(p) = gP 7 n 74i + -p^, 



(34) 

(35) 
(36) 
(37) 
(38) 
(39) 



where p 7 and p e (= m e n e ) are energy densities of photons and electrons, and v\ are their bulk three velocities 
defined by v % = u l /u°, and IT!^ is anisotropic stress of photons. It should be noted that the collision term ([34]) was 
obtained nonperturbatively with respect to density perturbations [30I ] . 



B. Evolution Equations of Magnetic Fields 



Now we obtain the Euler equations for protons and electrons as 

0, 

4orp 7 an 



TTlpTLUpUpi.^ 



enu£F ifl 



(Vei -V ji ) + jVej^i 



(40) 
(41) 
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where m p is the proton mass. Here we ignore the pressure of proton and electron fluids. Also the Coulomb collision 
term is neglected as explained below Eq. ([5]) . Note that the collision term was not evaluated in a manifestly covariant 
way. Here the left hand side in Eqs. ([40l) and (j4Tjl should be evaluated in conformal coordinate system. We also 
assumed the local charge neutrality: n = n e ~ n p . In the case without electromagnetic fields (i 7 ^ = 0), the sum of 
the equations (|40[) and (|4Tj) gives the Euler equation for the baryons in the standard perturbation theory. On the 
other hand, subtracting Eq. (|40p multiplied by m e from Eq. (|41|) multiplied by m p , we obtain 



m p m e 



m e en 



- Ui 



+ en{m p + m e )u^Fi^ - (m p - m e )j fl Fi f 



Ampp^anaT 



(42) 



where u M and j M are the center-of-mass 4- velocity of the proton and electron fluids and the net electric current, 
respectively, defined as 



m p u p L + m e u^ 



nip + m e 



en(u p 



(43) 
(44) 



Employing the Maxwell equations F£" = j^, we see that the quantities in the square bracket in the l.h.s. of Eq. (|42|) 
is suppressed at the recombination epoch, compared to the second term, by a factor 35] 



3 x 10 



_ 40 /10 3 cur 3 \ / IMpc 



(45) 



where c is the speed of light, L is a characteristic length of the system and oj p = \J Amne 1 1 m e is the plasma frequency. 

The third term in the l.h.s. of Eq. (|42|) . i.e., (m p — m e )j^Fi ill is the Hall term which can also be neglected because 
the Coulomb coupling between protons and electrons is so tight that \u l \ ^> \vL — u\\. Then we obtain a generalized 
Ohm's law: 



3e 



V 1% ) + -jVejlLj 



(46) 



Now we derive the evolution equation for the magnetic field, which can be obtained from the Bianchi identities 

F\p. v ,\] = 0, as 



= - e «* u /*F b . M 



2 

u^B 1 



3,,* 



B l v?), 



(47) 



where e y& is the Levi-Civita tensor and B l = (a 2 B l ) = e t: > k Fjk/2 is the magnetic field in the comoving frame [36j. We 
will now expand the photon energy density, fluid velocities and photon anisotropic stress with respect to the density 
perturbation as 



(") (i) 
p y (t,Xi) =Pt(t)+ Pj{t,Xi) + ■ ■ 

u l (t, Xi ) = u l (t,x l ) + { u l (t, Xl ) + 
T&{t, Xi ) =n«(t,Xi) + --- , 



u°(t,Xi) = a(t)~ 1 + u°(t,Xi) H 

Vi(t,Xi) = v t (t,Xi)+ Vi(t,Xi) + 



(48) 



where the superscripts (0), (1), and (2) denote the order of expansion and t is the cosmic time. Remembering that B l 
is a second-order quantity, we see that all terms involving B l in Eq. (|47|) . other than the first term, can be neglected. 
Thus we obtain 



d& 
~dt 



<0) 

4ctt P yd 



3e 



7" ijk 



( _W \ /(i) (i) \ l/(i) (L), \ „ AD (i) \ /(2) (2) \ 

^ 7>J — 2 [v ek - v lk j -- \^v el n 73 J -2 h j, k { v ei- - {v ejik - u 7J - fe J 



49) 
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where $and h l j are first-order curvature and tensor perturbations, respectively, in Poisson gauge, and we used the 

(!) (1) (0) 

density contrast of photons, 5 -y,k =P^,kl Py Further, we employed the fact that there is no vorticity in the linear 

order: e yfc v j /. — 0. It should be noted that the velocity of electron fluid can be approximated to the center-of-mass 

velocity at this order, v \ ~ v i . The physical meaning of this equation is that electrons gain (or lose) their momentum 
through scatterings due to the relative velocity to photons (baryon-photon slip), and the anisotropic pressure from 
photons. The momentum transfer from the photons ensures the velocity difference between electrons and protons, 
and thus eventually generates magnetic fields. We found that the contribution from the curvature perturbation is 
always much smaller than that from the density contrast of photons in the first term in Eq. (|49|) (which will be clearly 
seen in FigsfT] and Therefore we shall omit the curvature perturbation hereafter when considering the evolution 
of magnetic fields. 

The first term in Eq. (|49|) is exactly the same discussed in [29]. They have estimated contributions from these 
terms by considering typical values at recombination. Here we solve the equation numerically and obtain a robust 
prediction of the amplitude of magnetic fields. 

Eq. (|49]) shows that the magnetic field cannot be generated in the first (linear) order. The r.h.s. of Eq. (|49]) contains 
two kinds of source terms, i.e., intrinsic second order quantities and products of first order quantities. Since the first 
order quantities can be exactly evaluated within the frame work of the standard cosmological linear perturbation 
theory, we hereafter concentrate on the products of first order quantities. The evaluation of the intrinsic second order 
term will be left for the future work. 



III. COSMOLOGICAL MAGNETIC FIELDS 



A. Spectrum of Cosmologically Generated Magnetic Fields 

In this section we derive the spectrum of magnetic fields generated by cosmological perturbations. The cosmological 
perturbations can be decomposed into scalar, vector, and tensor modes. These modes physically correspond to 
perturbations of density, vorticity and gravitational waves, respectively. These cosmological perturbations are most 
likely to be generated during the inflation epoch. Most of the simple single field inflation models predict generation of 
adiabatic scalar and tensor type perturbations but vector type perturbations. Even if the vector type perturbations 
are generated, they are known to be damped away in the expanding universe. Observations of CMB temperature 
fluctuations and large scale structure of the universe strongly suggest the adiabatic scalar perturbations while we 
only have upper limits on tensor perturbations [37| . Therefore at most tensor perturbations are only sub-dominant 
components. Hence we consider only scalar type perturbations for linear order quantities throughout this paper. 

The linear-order scalar perturbations can be written as 

d 7 , k = 7=L= f(ik k )6 7 (k,ty ts d% (50) 



Vel = _L_ l{-ik)v & {k,ty^d% (5i) 



v el , k = -== I kk k v e (k,t)e i,! - s d% (52) 
V( 27r ) 



H< = — L= f(-k% + U)U(k,t)e^d 3 k, (53) 
v/(2tt) 3 J 3 

U l u = -j±= J{ik k ){-k% + \^n{k,ty^d 3 %. (54) 

Here k = k l jk. By using these expressions, we can write the evolution equation of magnetic fields in Fourier space 
as, 



(a*B\K,t)y = ^i£i e m—L= I d zg 

y y " 3e v/(2^f 



k'jKkS^K - k', t)5v(k', t) ~ -a ]k (K, k')v e (K - k',t)U(k', t) 



(55) 
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where we have defined, 



Sv 

a jk (K,k>) 



v e - v 
K 



7> 



\K-k'\ \ 

To obtain the spectrum, one needs to evaluate (B i (K)B* (K')) = S(K)S(K - K'). That is 



(56) 
(57) 



a 4 B l (K)B*(K>) = 



ijk Im 



(2 



^ [ d 3 k'd 3 k [° dt' [° dt" 

J Jo Jo 



k'.KkS^K - k' ,t')Sv(k' ,t')a(t') Pj (t')kiK' m 6*(K' - k,t")6v*(k,t")a(t") Pl (t") 
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a jk (k,k')v e (K - k\t')n{k\t')a{t') Pl (t')ai m (k',k)v* e {k' - k,t")Yl*(k,t")a(t") Pj {t") 



58) 



The final task we should do is to take an ensemble average of this expression. Assuming that perturbations are 
random Gaussian variables and using the Wick's theorem, we can expand the ensemble average of the products of 
four Gaussian variables as, for example, 



S^K - k', t')Sv{k', t')6*{K' - k, t")Sv*(k, t") 



S^K - k',t')6v{k',t')) U*jk> - k, t")6v*(k, t") 



+ (d^K - k>,t')5*(K> ~ k,t")j {5v{k',t')5v*{k,t") 
+ (s 7 (K - k',t')6v*(k,t")\ (Sv(k',t')51(k> - k,t") 



Since the evolution equations for linear order perturbation variables are independent of k, we may write 

5j(k,t) = tpi(k)S 7 (k,t) , 
Sv(k,t) — xpi{k)8v(k,t) , 
n(M) = ^(k)U(k,t) , 



(59) 

(60) 
(61) 
(62) 



where ipi{k) is the initial perturbation and 5 1 {k, t),Su(k, t), and II(fc, t) are the transfer functions, which are the solu- 
tion of the Einstein-Boltzmann equations under the adiabatic initial condition with ipi(k) = 1. By virtue of linearized 
equations, the overall amplitude and spectrum of the primordial fluctuations, P{k), can be given independently of 
the transfer functions. Here P(k) is defined using the two-point correlation function of tpi in Fourier space, 

1 p i (k)^*(k')) = P(k)5(k-k') , 



(63) 

where <5 is the Dirac delta function. According to the convention, we write P(k) oc /c" s ~ 4 with n s being the spectral 
index of primordial density perturbations. Using these definitions, we finally obtain, 



S(K) 
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V 3e J J (2tt) 3 
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'L T S 1 (k,\K-k\)S 1 (\K-k\,k) 
\K - k\ 
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/3S 2 (k, \K - k\)Si{\K - k\,k) 



— -aSAk, \K - k\)S 2 {\K -k\,k) 
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aSi(\K - k\, k) - pS 2 (k, \K - k\)S 2 (\K -k\,k) 



(64) 




FIG. 1: Evolution of perturbations relevant with magnetic field generation in a ACDM model with wavenumber k/h — f.6 
Mpc~ . Left: the five lines represent density contrast of CDM (<5cdm; black dotted line), and density contrast (<5 T ; red 
line) , velocity perturbation (t> 7 ; blue dash-dotted line), anisotropic stress (IT 7 ; green dashed line) of photons, and curvature 
perturbation ($; brown thin line) as indicated. Right: time evolutions of density contrast (8b\ red line) and velocity perturbation 
(vb; blue dash-dotted line) of baryons. The baryon-photon slip (v~, — v e ; blue dash-dotted line) is also shown in the figure. In 
the deep radiation dominated era, baryon-photon slip and anisotropic stress of photons are severely suppressed compared to 
the other perturbation variables. Perturbation variables are shown in the Newtonian gauge [3^ while anisotropic stress and 
slip terms are gauge invariant. 



where 

Si (MO = / ' a{t')p^{t')5^{k,t')5v{k',t')dt' , (65) 

5*2 (MO = / ' a(t')p y (t')v e (k,t')Il y (k' ',t')dt' , (66) 
3K(K -k)-2k 



(67) 



k {k -k{K -k)\ 2 

P = =, L - - ■ (68) 

\K-k\ 2 3 v ' 

Here Si and S2 describe effects of baryon-photon slip and photon anisotropic pressure. Specifically, the terms in the 
second and fourth lines in Eq. (|64|) are what we will call in subsequent sections the contributions to the generation of 
magnetic fields from baryon-photon slip and anisotropic stress of photons, respectively. 



B. Cosmological Perturbation Theory and Numerical Calculations 



Evolutions of perturbations relevant to the generation of magnetic fields can be solved within the frame work of 
the standard cosmological perturbation theory. The standard theory describing evolutions of linear perturbations, 
which was originally developed by Lifshiz [3^| and summarized in his textbook and other recent articles jU H2| , 
is well tested by a number of observations and firmly established. As for an application to cosmological models, the 
first realistic numerical calculation was done by Peebles & Yu [43j, where the distribution of photons was solve d by 
directory integrating the Boltzmann equation. The calculations were subsequently extended to include neutrinos [4J], 




k/h [Mpc -1 ] 



FIG. 2: Spectra of <5cdm (black dotted line), <5(, (thick black line), <5 7 (thin red line), II 7 (green dashed line), Vb — u 7 (blue 
dash-dotted line), and <& (brown dash-dot-dotted line) at redshift z = 10 3 (bottom left), 7 x 10 3 (top right), 5 x 10 4 (top left). 
In general, the spectra at each redshift are divided into four parts from large scales to small ones: primordial spectra at super- 
horizon scales, acoustic oscillation spectra at sub-horizon scales, exponentially damping spectra at diffusion scales, and power 
law phase after the diffusion damping. In acoustic oscillation scales at each redshift, the amplitudes of density perturbations 
of baryons (Si,) and photons (<5 7 ) stay constant, while baryon-photon slip (v e — « 7 ) and anisotropic stress of photons (II 7 ) show 
larger power at relatively smaller scales. Powers in baryon-photon slip and anisotropic stress of photons get larger at lower 
redshift. 



and non-baryonic dark matter [45| . These calculations, especially the resulting angular spectrum of CMB photons, 
are fully understood within an analytic treatment by Hu & Sugiyama (4o| . 

Evolutions of the linear quantities such as 5, v, and II in an exp anding universe can now be easily solved by the 
publicly available Einstein-Boltzmann code such as CMBFAST j47j or CAMB (48|. In Fig. [TJ we show an example of 
the evolution of perturbations with wavenumber k/h = 1.6 Mpc -1 in standard ACDM cosmology. Time evolutions 
can be typically divided into four characteristic eras. First, the modes of perturbations are beyond the cosmic horizon 
(era I; a < 10~ 5 in Fig. [J). After the horizon crossing, the baryon-photon fluid undergoes acoustic oscillations (era 
II; 10~ 5 < a < 10~ 3 - 5 ) until the diffusion of photons erases the perturbations (era III; 10~ 3 5 < a < 10~ 3 ) as the 
universe expands, while density perturbations of CDM grow mildly by their own gravity. After recombination (era 
IV; a > 10 -3 ), photons and baryons are decoupled each other. In that era photons stream freely and baryons can 
evolve according to the gravitational potential of CDM. 

All modes of perturbations in baryons and photons whose wave number is larger than ~ 0.1 Mpc -1 arc destined 
to be erased before recombination by Silk damping [49j | . Therefore, each perturbation mode with such wavenumber 
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FIG. 3: Plots of integrand per dlog(l + z), k 6 ^5~,{u, 



(left) and k 4?IXyM e (right). The magnetic field spectrum results 



from non-linear convolution and time integration of these spectra. The contributions mainly comes from the waves between 
the scales of cosmic horizon, and the Silk diffusion at each redshift. 
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FIG. 4: Time evolution of Spectra of conformal magnetic fields generated from baryon-photon slip (left) and anisotropic stress 
of photons (right). Corresponding redshifts are also indicated in the figures. For both contributions, the smaller scale fields are 
created at the higher redshift. 



can contribute to the generation of magnetic fields only when the mode undergoes acoustic oscillations. During that 
era, density and velocity perturbations of baryons and photons show an oscillatory behavior and their amplitudes 
remain almost constant. Meanwhile the variables suppressed by tight coupling between photons and electrons, such 
as v e — or Il 7 , grow as the universe expands and the number density of electrons becomes smaller (see Fig. [TJ. 
Thus the spectra of perturbations have a peak at the scale of acoustic oscillations or the scale slightly larger than the 
diffusion scale at each epoch, that is clearly shown in Fig. O 

We plot {ap 1 /H)k :i P{k)5v{k)8 1 {k) and {ap 1 /H) k 3 P(k)v b (k)U 7 (k) in the left and right panels of Fig. H These 
spectra naively correspond to the source terms S\ and S2 in Eq. (|64[) . Of course, the magnetic field spectrum should 
be obtained by a non-linear convolution of these spectra given by Eq. (I64p . However, because the smaller scale modes 
than the diffusion scale are absent at each redshift and thus the non-linear couplings mainly come from the modes 
with comparable wavenumbers, one can expect that the cosmologically generated magnetic fields will roughly have a 
spectrum similar to the envelope curves of these spectra. 

We numerically integrated Eq. (|64j) to obtain the magnetic field spectra. The cosmological parameters in our 
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calculation are fixed to the standard A CDM values in a flat universe, i.e., 

(h,n s ,Q b ,n A ,A) = (0.7, 1.0, 0.04, 0.70, 2.7 x 10~ 9 ) , (69) 

where h is the current Hubble parameter in units of lOOkm/s/Mpc, n s is the power spectrum index of primordial 
density fluctuations, flf, and f^A are energy densities in baryon and cosmological constant in units of the critical 
density, and A is overall amplitude of density perturbation squared. Magnetic field spectra induced by baryon-photon 
slip and anisotropic stress of photons at different redshifts, which correspond to second and fourth lines of Eq. ([M|) . 
respectively, are shown in left and right panels of Fig. [U In these figures we take the combination of y^k 3 (BB*) /2n 2 
to show in units of Gauss. The spectra consist of two slopes. At super-horizon scales the slope is proportional to fc 3 5 
for both contributions while it becomes less steep at sub-horizon scales. The turning points of the spectra roughly 
correspond to the Hubble horizon scale at each redshift while it may be difficult to see clear turnoff in the spectra 
generated from anisotropic stress (right panel). Clearly, one can find that small scale magnetic fields are created at 
higher redshifts. After their generation, magnetic fields B are adiabatically diminishing as a~ 2 . To see this fact more 
explicitly, we depict the time evolution of magnetic field with fixed wavenumber K in Fig. [5l Analytic interpretation 
for these magnetic field spectra is given in the next section. 



IV. ANALYTIC INTERPRETATION OF THE MAGNETIC FIELD SPECTRUM 



A. Magnetic Field Spectrum at Large Scales 



Magnetic fields are mainly created after the modes of perturbations with the corresponding scale enter the cosmic 
horizon and are causally contacted. Magnetic fields at larger scales than the cosmic horizon are generated only as a 
consequence of non-linear couplings of each Fourier mode. In the previous studies it was suggested that at largest 
scales the magnetic field spectrum generated from density perturbations has the power proportional to k 2 [27| ■ If one 
carefully evaluates, however, the terms proportional to k 2 exactly cancel each other out, and the spectrum starts from 
the term proportional to fc 4 as we shall show below. In order to see this, let us consider the spectrum generated by 
the first term (baryon-photon slip term) in equation (|49p . 



B l (K,t)B*(K',t)^ oc J d 3 kK 2 {l - (jfc ■ Kf] P{\K - k\)P{k) 



^SfQk-kl^-j^^^ikAk-k^Qk-klk)} . (70) 



In the limit K/k — > 0, we can write 



\K - k\ = y'k 2 - 2kKj + K 2 fak-Kj , 

\k-k\- 1 « k~ l (i + ji) , 

where 7 = k • K. Furthermore, time integration terms can be treated as K independent in the same limit, i.e., 
S\(\K — k\,k) w S\(k, \K — k\) ps T(k). Remembering that the spectrum of density perturbations is written as 
P(k) cx fc" s ~ 4 , we can rewrite equation ([70)) as 

B l {K,t)B*{k',t)^ cx J k 2 dk J d 7 (l - j 2 )P{\K - k\)P{k)K 2 ^1 - ^ T 2 {k) 

,2j7. / a\ r , / aw^a 1.71,-4 / n t2/ 



k 2 dk J d 7 ( 7 - r ) [k - (n s - 4)Ki\ k n ^ 4 I — —\ T 2 (k) 

cx K 4 . (71) 
We can see this non-linear power law tail at large scales in Fig. [U 



B. Magnetic Field Spectrum at Sub- horizon Scales 



In this subsection we show that the behavior of the magnetic field spectrum at small scales can be understood within 
the standard theory of cosmological perturbations. Since Compton scattering is an essential process for generation 



scale factor 



scale factor 



FIG. 5: Time evolution of magnetic fields in units of Gauss, ^Jk 1, S(k) /2ty 2 , generated from baryon-photon slip (left) and 
anisotropic stress of photons (right) with fixed wavenumbers as indicated. Magnetic fields at smaller scales have been created 
earlier and their amplitudes are larger. After their generation magnetic fields diminish as a -2 adiabatically. 



of magnetic fields, magnetic fields start to be generated when the relevant modes come across the horizon and are 
causally connected. On the other hand, once the modes become shorter than the diffusion scale of photons, magnetic 
fields cannot be no longer generated. Therefore magnetic fields with wavenumber K are mainly created from density 
perturbations with wavenumber k ~ K < fedis, where k^ie is the corresponding wavenumber of the diffusion scale at 
each time. Within the approximation that k ~ K, \K — k\ ~ K, the spectrum of magnetic fields is roughly given by 




where 



si = P{k)S 1 {k,k) , (73) 
s 2 = P(k)S 2 (k,k) . (74) 

Here we neglected sub-dominant cross correlating terms. We have already plotted k 3 s± cx k 3 ^j-P(k)5v(k)8 1 and 
k 3 S2 <x k 3 P (k)vb(k)H 7 in the left and right panels of Fig. [3) These spectra indeed indicate that magnetic fields 
are created when the modes of perturbations come across the cosmic horizon and undergo acoustic oscillations. In 
fact, the largest contributions come from the modes of perturbations slightly larger than the diffusion scale. 

The behavior of acoustic oscillations in the early universe can be analytically understood by using tight coupling 
approximations. In the tight coupling regime, Compton scattering occurs sufficient enough so that photons and 
baryons (electrons) can behave as a coupled fluid. Therefore, the baryon-photon slip (vb~ i> 7 ) and the anisotropic stress 
of photons II 7 are severely suppressed. Specifically, within the tight coupling approximation [46| the perturbation 
variables satisfy, 



4p 7 



where e is a tight coupling parameter roughly given as 



(75) 
(76) 
(77) 
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This parameter represents how tightly photons are coupled with electrons and gives a criteria whether one can treat 
photons and electrons as a single perfect fluid. So, the parameter in equations ([76]) and (177)) simply reflects the fact 
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that when the photons and electrons can be treated as a tightly coupled single fluid (i.e., when the tight coupling 
parameter e is sufficiently small), the baryon-photon slip (v e — i> 7 ), and the anisotropy of the distribution function 
(IXy) are severely suppressed in comparison with the density fluctuation (<5 7 ) (see Fig. [T]). 

Let us see how the spectra in Fig. [3] can be understood using the above relations. Noting that Hubble parameter 
scales as H oc a~ 2 in the radiation dominated era, we have the following relations 

JjriVb-V^ OC -—^ ev l [j^J X kS l a ( k ) > ( 79 ) 

— 4 

-jll 7 oc — ^-cu-v oc ka°S 1 , (80) 
H a 

where we have used the tight coupling relations (|75"|) , (|75|) . and (|77|) . In the case where the spectrum of primordial 
density perturbations is scale invariant, k 3 P(k)65 is constant when the mode with wave number k enters the Hubble 
horizon. Therefore, the spectra in Fig. [3] should be proportional to 

fc 3 si oc k 3 ^P(k)S 7 (v b {k) ~ v 7 (k)) oc ka 2 , (81) 
H 

k 3 s 2 oc /c 3 P(fc)n 7 i; b ^ oc ka . (82) 
H 

These behaviors are clearly seen in Fig. [3J 

However, a complication arises when considering the spectrum of magnetic fields. Magnetic fields can not be gen- 
erated from the first order solution of tight coupling approximation [501 ]. For example, one finds that the contribution 
from the cross product between the gradient of density perturbation of photons (<5 7J ) and velocity differences of elec- 
trons and photons (v e j — u 7 j) (the first term in Eq. (l49p ) vanishes because these vectors are colinear in the first order 
tight coupling approximation. Therefore, we should consider the second order solutions of tight coupling expansion 
which are proportional to 0(e 2 ). Accordingly, the part of the spectra which can contribute to magnetic fields (Eqs. 
(fBTj) and (|52")l ) should now be evaluated as 

fc 3 si|ma g oc k 3 '^-P{k)5 1 {v b {k) - u T (fc))| mag oc e x fca 2 oc k 2 a 4 , (83) 
H 

fc 3 s 2 |ma g oc fc 3 P(fc)n 7 u b ^| mag oc e x ka oc k 2 a 3 . (84) 

From these relations we found that the contribution to magnetic fields becomes larger as the universe expands larger 
(the larger a) . This means that the scale factor a should be evaluated at the onset of the photon diffusion where the 
density perturbations are erased, because the contribution becomes largest at that time. The diffusion length is given 

by m 

fc-2. = I [ 1 i+ iffl ^ (85) 

dlff 6 J an e a T 1 + R a ' V ' 

from which one obtains the relation between the diffusion length and scale factor in the radiation dominated era: 

fcdiff = 4.3 x (j^r) 3 2 Mpc~ x oc a~ 3 / 2 . Combining these considerations altogether, the spectrum of magnetic fields 
at sub horizon scales is roughly given by, 

S(k) - k 5 sl + k 5 s 2 2 

- £T 7/3 (slip term) + AT 1 (II 7 term) . (86) 

At smaller scales anisotropic stress of photons (n 7 ) should be the dominant contributor to the magnetic fields. 



C. Cut-off of the magnetic field spectrum 



As we discussed in the previous subsection, the strength of magnetic fields increases as B ~ ^k 3 S(k) oc k on small 
scales where the contribution from anisotropic stress of photons is dominant. One might expect there exists more and 
more magnetic fields on smaller scales. However careful treatment is needed for the generation of magnetic fields on 
very small scales since the relativistic effects of electrons, which we omit in the formulation developed in section II, 
play an important role when the temperature of the universe was higher than the electron rest mass, i.e., T > 511keV. 
The power spectrum of magnetic fields below the diffusion scale at T — 511keV, therefore, would be modified. 
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There are two potential effects which affect the generation of magnetic fields quantitatively when the temperature 
of the universe was higher than the electron rest mass and electrons were relativistic. First one is the transition 
of scatterings between photons and electrons from the Thomson regime to the Compton one. The cross section 
of Compton scattering is proportional to the inverse square of the center of mass energy while that of Thomson 
scattering is independent of the energy. At first glance, this transition would result in the reduction of the electric 
current and hence the reduction of magnetic fields because photons can not push electrons effectively (Eq. |49|) . 
Interestingly, however, the weaker interaction makes the velocity difference between photons and electrons and the 
anisotropic stress of photons larger, and these two effects would cancel each other out. Specifically baryon-photon 
slip and anisotropic stress of photons responsible for the generation of magnetic fields are proportional to the tight 
coupling parameter e in equation (|78[) . which contains <7t in the denominator. Inserting this into equation (|49[) . we 
find that magnetic fields generated from cosmological perturbations would be cross section independent, i.e., 

B oc a% . (87) 

The second relativistic effect to be considered is the electron-positron pair creation from photons through 27 <-> 
e + +e~. In the thermal history of the standard big bang paradigm, electrons and positrons were in thermal equilibrium 
with photons before the temperature of the universe was around and higher than the electron rest mass. At that time 
their number densities were about the same as that of photons. As the universe expanded and cooled, photons became 
less energetic and unable to pair-create electron-positron pairs while existing electron-positron pairs annihilated into 
photons. The number density of electrons decreased drastically when the universe cooled below the critical temperature 
corresponding to the electron rest mass, and consequently, tiny residual number of electrons compared to that of 
photons, n e /n 7 ~ 10~ 9 , has survived. In other words, if one looks back in the thermal history of the universe, 
the number density of electrons increased all of sudden by a factor of 10 9 around the critical temperature. This 
effect would significantly modify the magnetic field spectrum. Since couplings between photons and electrons (and 
positrons) were so tight above the critical temperature than below, the relative velocity between them or anisotropic 
stress of photons, and thus induced electric current, were smaller above the critical temperature. Specifically, the 
tight coupling parameter is proportional to the inverse of electron number density, so are the magnetic fields, 

Ba^ 1 , (88) 

which leads to the reduction of the amplitude of the magnetic field spectrum by a factor of 10 -9 at the corresponding 
scale. 

The diffusion scale at the critical temperature (T — 511keV) is ~ 2 x 10~ 3 pc in the comoving scale. Therefore we 
expect that the spectrum of magnetic fields continues to increase up to this scale and rapidly drops on this scale. There 
exists negligible amount of magnetic fields on smaller scales. To summarize, the magnetic field spectrum generated 
from scale-invariant density perturbations is given as 

fc 3 5 ( A > 100 Mpc) 

B oc oc I fcV3 (1 ° MpC < A ~ 100 Mpc) (89) 

i * 0C V**W«S k (io-3 pc <A<10 Mpc) 

- ( A < 10~ 3 pc) 

We depict the spectrum of magnetic fields in Fig. [51 By remembering the fact that the spectrum of magnetic field 
S(k) is proportional to P(k)P(k — fc'), the spectrum can be easily generalized in the case of tilted spectrum: 

{k 3 5 ( A > 100 Mpc) 

fc-" 2 / 3 (10 Mpc<A<100Mpc) 

k n ° (10- 3 P c<A<10 Mpc) (9Uj 

- ( A < 10~ 3 pc) 

where n s is the spectral index of primordial density perturbations. 



V. MAGNETIC HELICITY 



So far, we have investigated the spectrum of magnetic field amplitude while there exists another important quantity, 
helicity, which defines the property of stochastic magnetic fields. Magnetic helicity is a quantity which characterizes 
magnetic field configuration how much the magnetic field lines are twisted, or how many closed magnetic lines are 
linked [51| . Among the cosmological mechanisms to generate magnetic fields, a number of scenarios predicts the 
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FIG. 6: Spectrum of magnetic fields generated from cosmological density perturbations at z = 1100. Red dashed line shows 
contribution from baryon photon slip, blue dash-dotted line from anisotropic stress of photons, and black like shows the total 
spectrum. The spectrum was obtained by direct numerical integration up to k=30 Mpc -1 , and by analytic extrapolation above 
it. The spectrum has a slope proportional to fe 3 ' 5 at larger scale than the horizon at recombination. At smaller scales than 
~1 Mpc -1 , magnetic fields are dominated by the contribution from the anisotropic stress of photons. At smallest scale around 
2 x 10 -3 pc, the spectrum would have a cut-off due to the relativistic effects of electrons (and positrons). 

primordial fields with non-zero helicity [52l . l53l . [54|. Therefore, helicity, as well as amplitude, can be an important 
probe of the origin of primordial magnetic fields. 

Since magnetic helicity is a conserved quantity in the limit of infinite conductivity in the standard MHD theory, it 
should remain zero if there was no helicity initially in the early universe. To generate cosmological magnetic fields with 
non-zero helicity, parity violating processes would be necessary, often related with CP violation of fundamental particle 
interactions. This will lead us to the conclusion that helicity can not be generated if magnetic fields were generated 
through density perturbations where the standard Compton scatterings are the only relevant particle interactions. 
In fact, we can explicitly show that magnetic fields generated from cosmological perturbation discussed in this paper 
have no magnetic helicity as follows. 

Formally, the magnetic helicity density is defined as 

H = y J cPxA-B , (91) 

where A is the vector potential and V is the normalization volume. Similar to the way to define the magnetic field 
spectrum S(k), the helicity spectrum Pn(k) can be defined as 

(Mk)B*(k>)) = 5^{k-k')P ll (k)^> , (92) 

where Pij(k) = Sij — kikj is the projection tensor Note that Py satisfies the divergence-less condition of magnetic 
fields, k*Pij = 0. Remembering that Bi = ieijkkjAk, Eq. (|92p can be put together with magnetic field spectrum S(k) 



[f] In general, Pij is not necessarily a symmetric tensor. In fact, three conditions we put on magnetic fields, i.e., reality, divergence-less, and 
statistical homogeneity and isotropy require that Pij should have the form: Pij{k) = f(k) u5y — fe'fc'l + ig(k)e i jik l l55l l56( . However, 
as long as magnetic helicity is concerned, antisymmetric part in Eq. 1921 1 is irrelevant because helicity is defined as an inner product. 
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as 



B i (k)B*(k'))=6^(k 



P tJ {k) 



S{k) 



+ itijik 



P H {k) 



(93) 



Thus, the antisymmetric part of magnetic field correlations denotes helicity of magnetic fields. 

To see explicitly that magnetic fields generated from density perturbations have no helicity, let us consider an 
alternative expression for Eq. 



Bj(k)B* (fc')) = 5^{k-k') 



P tJ (k) 



S(k) 



p H (k) 



(94) 



and subtract this from Eq. (|93p , which should be zero if there is no magnetic helicity. The magnetic fields generated 
from cosmological perturbations can be symbolically expressed as 



Bi(k) 



Aab 



d 3 KK a k b f{K, k)S(k - K)^(K) 



(95) 



where f(K, k) is a real function of K and k (eq. (|55[) ). and 8 and 7 are the time integration of perturbation variables. 
Let us now explicitly evaluate I^Bi{k)B*(k')^ — (^Bj(k)B* (fc')^ as 

Bi{k)B*{k')\ -/Bj(k)B*(k')) = 6 (3) {k-k') f d 3 KP(K)P(\k - K\)f(K,k)e lab e : > lm 



x {K a k b K t k m - Kik m K a k b )f(K, fc) 7 2 (X)<5 2 (|fe - K\) 



+ 



(k a k b (ki - Ki)k m - Kik rn (k a - K a )k b 



\k-K\ 



= . 



f(k ~ K, k)S(\k - K\)j{\k - K\)j(K)5(K) 

(96) 



Thus, the helical part of magnetic fields from density perturbations should vanish. This fact is in marked contrast 
to the other origins of primordial magnetic fields in the very early universe (above GeV), where the magnetic fields 
must be somewhat helical due to the interaction of the magnetic field with a cosmic axion field • 

Several ways to detect helicity of magnetic fields have been proposed. It was shown that primordial helicity can be 
a new source to induce parity-odd correlations such as between temperature anisotropy and B-mode polarization, and 
that between E-mode and B-mode polarizations in CMB anisotropies [5 a , l58l . |59| , which are zero for fields without any 
helicity. Recently, Kahniashvili and Vachaspati proposed that the correlation of the arrival momenta of the cosmic 
rays can be used to detect helicity of an intervening magnetic field [6(| ■ Although it is very difficult to detect helicity 
in the cosmological magnetic fields, yet it deserves further investigations since detecting helicity will shed light on the 
mystery of the origin of large-scale magnetic fields. 



VI. DISCUSSION AND SUMMARY 



In this paper we discussed a generation mechanism of magnetic fields from cosmological density perturbations. 
Following our previous papers, we present a detailed formulation for the magnetic field generation. The key is that 
photons scatter off only electrons (and not protons) by Compton scatterings. Electric fields are induced to prevent 
charge separation between electrons and protons, and magnetic fields are generated from these induced electric fields 
through Maxwell equations. We showed that the baryon-photon slip and the anisotropic stress of photons can generate 
the magnetic fields if second order couplings in density perturbations are taken into account when evaluating Compton 
scattering terms. Since these two sources for magnetic fields naturally arise from density perturbations for a wide 
range of wavenumbers, the fields are also naturally generated for a wide range of scales. We then gave an analytic 
interpretation of the resultant spectrum of magnetic fields within the framework of the cosmological perturbation 
th eory. U sing the tight-coupling approximation of baryon-photon plasma, we showed that the fields have the spectrum 
\J k 3 S(k) oc k at small scales where the photon anisotropic stress is the dominant contributor, and thus the fields 
become stronger at smaller scales. A typical amplitude of magnetic fields is of the order 10 -24 Gauss at 1 Mpc scale. 
Magnetic helicity should not be associated with these magnetic fields. We also found that magnetic fields have a 
cut-off around 2 x 10~ 3 pc due to the relativistic effects of electrons. 
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In our previous papers [3(| [Ml we suggested that the magnetic fields as strong as 10 -18 Gauss can arise at Mpc 
scale at recombination, which is larger than the amplitude presented in this paper. The main difference from the 
previous results comes from the missing scale factor in equation (|49|) . This mistake leads an overall overestimate for 
the magnetic field amplitude and a steeper magnetic field spectrum on small scales, where the fields are mainly created 
at earlier epochs. The amplitude obtained here is widely consistent with recent literatures around a comoving scale of 
0.1 Mpc -1 [I?], [H, [6l[ . However, there still exists the difference of the amplitude of magnetic fields at smaller scales. 
For example, our results show stronger magnetic fields by the factor of 10 2 than those reported in [U [6l| at Mpc 
comoving scale, and the difference becomes even larger at smaller scales. This may be because they have evaluated 
the strength of magnetic fields at the instance of recombination at which the Silk damping effect have already erased 
perturbations on scales smaller than k ~ 0.1 Mpc -1 , while we have solved the evolution equation of magnetic fields 
from deep in the radiation epoch, through matter-radiation equality, to recombination. In ref . [H[ , it is also discussed 
that larger magnetic fields arise at small scales if one consider the earlier period. It is natural to consider that magnetic 
fields generated prior to recombination can survive, because the diffusion scale in the highly conductive primordial 
plasma due to the Ohmic dissipation is much smaller than cosmological scales considered here. Therefore, we think 
that magnetic fields generated from (nearly scale invariant) density perturbations should dominate on small scales. 
We will give further details about the discrepancy between our work and others elsewhere soon [501 ] . 

As discussed in this paper, magnetic fields have a small scale power up to 2 x 10~ 3 pc comoving scale. This scale 
corresponds to the cosmic horizon when the temperature of the universe was around MeV. Until this epoch, thermally 
created electrons and positrons significantly suppressed the sources of magnetic fields. If we assume the scale invariant 
spectrum for primordial density perturbations, our results indicate that magnetic fields as strong as ~ 10~ 15 Gauss 
can arise at this scale. Although these are small compared with the micro Gauss fields observed in present galaxies, it 
would be enough for the hydro-dynamical dynamo action to amplify into the present amplitude during the structure 
formation until today Q. Of course, if density perturbations have a blue spectrum, or they are anomalously large at 
small scales [62J larger magnetic fields will be generated. 

Although cosmologically generated magnetic fields appear sufficient for 'seed fields' of galactic magnetic fields, it 
is still unclear whether they can act as seed magnetic fields for galaxy clusters. In most astrophysical objects, such 
as disk galaxies and stars, differential rotation is important for amplifying and sustaining their own magnetic fields. 
However, galaxy clusters have little rotation, and therefore, another mechanisms will be necessary. It is recently 
argued that cluster plasmas threaded by weak ma gne tic fields are subject to very fast growing instabilities, and this 
instability happens once the ions are magnetized [63:]. The magnetized condition corresponds to the magnetic field 
amplitude of 10~ 18 G for typical cluster parameters, which is a little larger than that of cosmologically generated 
seed fields which we found here. Furthermore, the fields amplified by plasma instabilities have a rather small reversal 
scale (typically 10 4 km to lOpc [64[), while the observational data suggest that the typical reversal scale is ~ 1 
kpc. Therefore we can not conclude that magnetic fields in clusters of galaxies can directly arise from the seed fields 
generated from density perturbations in the early universe. 

To observe these seed magnetic fields in a direct manner is interesting but very challenging. The magnetic fields at 
recombination should leave an imprint on CMB photons through Faraday rotation. Faraday rotation can then create 
B-mode polarization from the dominant E-mode polarization which associates with (scalar) density perturbations. 
Furthermore, sufficiently large magnetic fields can source the vector mode perturbations on baryon fluid by themselves, 
which then induce B-mode polarization pattern in CMB (65l . |66| . By analyzing polarization patterns at large scales 
one can in principle get rid of information about the magnetic fields. Observation forecast predicts that the fields as 
week as 10~ 10 G will be detected by future missions such as Planck, which is however still out of reach of the weak 
magnetic 'seed' field created from density perturbations discussed here. 

Another promising way to detect the weakest 'seed' magnetic fields would be to measure the time-delayed emissions 
from the gamma-ray bursts. The idea was originally proposed by Plaga (67j . and there the author concluded that 
magnetic fields as week as 10 -24 Gauss in void regions can be probed by detecting time-delayed events of GeV-TeV 
photons. Of course one may claim that it is highly uncertain whether magnetic fields in void regions purely consist of 
the remnant seed fields created in the early universe. Interestingly, it is argued that large volume of void regions can 
escape from any astrophysical activities such as galactic winds [681 ]. Therefore, it is likely that one can probe primordial 
magnetic fields by measuring void fields through gamma-ray photons. Unfortunately, it would be impossible to reach 
the field as weak as 10 -24 Gauss in practice, because there exist another effects to cause time delays to high energy 
gamma-ray photons. Subsequent studies have shown that magnetic fields strength as strong as 10 -19 G would be 
necessary for time-delays by magnetic deflection to be longer than those caused by angular spreading . Therefore, 
it appears that an amplification of 'seed' magnetic fields at recombination by at least ~ 100 will be necessary for the 
fields to lead any observational signatures in gamma-ray photons. Clearly it will be of great interest to study the 
evolution of magnetic fields in the void regions from recombination to the present universe. 
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